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We have developed a material specific theoretical framework for modelling scanning tunneling 
spectroscopy (STS) of high temperature superconducting materials in the normal as well as the 
superconducting state. Results for Bi2Sr2CaCu20s+5 (Bi2212) show clearly that the tunneling 
matrix element strongly modifies the STS spectrum from the LDOS of the dx2_y2 orbital of Cu. 
The dominant tunneling channel to the surface Bi involves the dx2_y2 orbitals of the four neigh- 
bouring Cu atoms. In accord with experimental observations, the computed spectrum displays a 
remarkable asymmetry between the processes of electron injection and extraction, which arises from 
contributions of Cu (i^2 and other orbitals to the tunneUng current. 

PACS numbers: 68.37.Ef 71.20.-b 74.50.+r 74.72.-h 



Scanning tunneling spectroscopy (STS) has entered the 
realm of high-temperature superconductors powerfully 
by offering atomic scale spatial resolution in combina- 
tion with high energy resolution. The physics of these 
materials is dominated by the cuprate layers, which are 
usually not exposed to the tip of the apparatus. Much of 
the existing interpretation of the spectra is based how- 
ever on the assumption that the STS spectrum is directly 
proportional to the local density of states (LDOS) of the 
Cu02 layer, neglecting the effects of the tunneling ma- 
trix element in the presence of the insulating over layers. 
Here, we focus on the Bi2212 system, which has been 
the subject of an overwhelming amount of experimen- 
tal workP, [2I, 0, 0, although our results bear more 
generally on the STS spectra of the cuprates. We take 
into account the fact that the current originating in the 
Cu02 layers reaches the tip after being 'filtered' through 
the over layers of SrO and BiO, and show that instead 
of being a simple reflection of LDOS of the Cu02 lay- 
ers, the STS signal represents a complex mapping of the 
electronic structure of the system. 

In order to construct a realistic framework capable of 
describing the tunneling spectrum of the normal as well 
as the superconducting state of the cuprates, we start 
with the normal state Hamiltonian 



Hi 



(1) 



which describes a system of tight-binding orbitals created 
(or annihilated) via the real-space operators cj^^ (or Caa)- 
Here a is a composite index denoting both the type of 
orbital (e.g. Cu (1^2 _y2) and the site on which this or- 
bital is placed, and a is the spin index, is the on-site 
energy of the a^^ orbital, a and P orbitals interact with 
each other through the potential Va/3 to create the energy 
eigenstates of the entire system. 

Superconductivity is included by adding a pairing in- 
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FIG. 1: (a) Side view showing tip placed schematically on 
top of seven layers used to compute the tunneling spectrum 
of Bi2212, where the surface terminates in the BiO layer. Tun- 
neling signal from the conducting Cu02 layers reaches the tip 
after passing through the filtering layers of SrO and BiO. (b) 
Top view of the surface showing the arrangement of various 
atoms. Eight two-dimensional real space primitive unit cells 
used in the present computations are marked by dashed lines. 



teraction term A in the Hamiltonian of Eq. [T] as follows 



H = Hi 



af3a 



(2) 



We take A to be non-zero only between (1^2 _y2 orbitals 
of the nearest neighbor Cu atoms, and to possess a 
d-wave form, i.e. A is given in momentum space by 
A/e = y [cos k^a — cos kya] , where a is the in-plane lat- 
tice constant. This interaction allows electrons of op- 
posite spin to combine into superconducting pairs such 
that the resulting superconducting gap is zero along the 
nodal directions kx = ^ky^ and is maximum along the 
antinodal directions. 

The Bi2212 sample is modeled as a slab of seven lay- 
ers in which the topmost layer is BiO, followed by layers 
of SrO, Cu02, Ca, Cu02, SrO, and BiO, as shown in 
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FIG. 2: (a) Normal state electronic spectrum of Bi2212. Cu d^2_y2 weight of states is given by the color of the dots. Inset 
shows the energy region of the VHS's on an expanded scale, (b) Electronic spectrum of the superconducting state. Electron 
weight of states is given by the color of the dots. Inset highlights the region of the superconducting gap near the antinodal 
point. 



Fig. 1(a). The tunneling computations are based on a 
2V^ X 2V^ real space supercell consisting of 8 primitive 
surface cells with a total of 120 atoms (see Fig. 1(b)). The 
coordinates are taken from the tetragonal crystal struc- 
ture of Ref. 0. A tetrahedral cluster of atoms attached 
to a thin slab simulates the tip and tip holder. The tip 
is allowed to scan across the substrate to generate the 
topographic STM map, or held fixed on top of a surface 
Bi atom for the STS spectra. 

The tight-binding parameters are fitted to the LDA 
band structure of Bi2212 that underlies for example 
the extensive angle-resolved photointensitycomputations 
of Ref. 7. The Slater-Koster results [sl, \^ are used to 
fix the angular dependence of the tight binding over- 
lap integrals. The specific orbital sets used for vari- 
ous atoms are: {s^Px^Py^Pz) for Bi and O; s for Sr; 
and (45, d^z'^_r-2, dxy, dxz^ dyz, dx2_y2) for Cu atoms. This 
yields 58 orbitals in a primitive cell, used in band calcu- 
lations, and a total of 464 orbitals for Green function 
supercell calculations in 256 evenly distributed k-points. 
Finally a, gap parameter value of |A| = O.OAbeV is cho- 
sen to model a typical experimental spectrum [1] for the 
generic purposes of this study. 

The LDOS and tunneling computations are based on 
Green function formalism. At first, the normal-state 
Green function is constructed via Dyson's equation using 
methodology described in Ref. Il0|. At this stage a self- 
energy for orbital a, = T^'^ i ^^^'a is embedded in 
Dyson's equation for possible effects of various bosonic 
couplings and correlation effects 11, For simplicity, 
we have assumed the self-energy to be diagonal in the 



chosen basis. In building up the Green function in the 
superconducting state, we utilize the conventional BCS- 
type self-energy ^^^^ = AG^At (see, e.g., Ref. [H), 
where is the hole part of normal state Green func- 
tion. 

Fig. 2 shows the calculated band structure of Bi2212 
in the normal as well as the superconducting state from 
Hamiltonians of Eqs. [T] and [2l The normal state is seen 
to properly display the major features such as: The pair 
of Cu02 bands crossing the Fermi energy (Ef) with the 
associated van Hove singularities (VHS's) marked VHS- 
a (antibonding) and VHS-b (bonding), split by 250 meV 
at the (tt, 0) point; BiO bands lying about 1 eV above 
Ef] and the 'spaghetti' of bands involving various Cu 
and O orbitals starting at a binding energy of around 1 
eV below Ef- Although states near Ef are mainly of 
Cu dx2_y2 and O Px,y character, they also contain some 
Bi and Cu d^^ admixture. In the superconducting state in 
Fig. 2(b), a quasiparticle spectrum mirrored through Ef 
is obtained with a doubled number of bands due to the 
pairing interaction. A d-wave superconducting gap opens 
up in both Cu02 bands near Ef- The quasiparticles have 
a mixed electron/hole character only near the edges of the 
gap. 

To compute the tunneling S2ectra we apply the 
Todorov-Pendry expression [ij, ll5| for the differential 
conductance a between orbitals of the tip {t^t') and the 
sample (5, s')^ which in our case yields 

= ^ = ¥ ^ Ptt'iEF)Vt,sP.AEF + eV)VJ,„ (3) 
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FIG. 3: (a) Computed (black line) and experimental (red line) differential tunneling spectra, dl/dV, as a function of bias 
voltage (in volts). Local density of states (LDOS) of Cu d^2_y2 electrons is shown (green line) as a function of energy, where 
the horizontal scale is in eV. (b) Computed topographic image of the BiO surface in which Bi atoms are bright and O atoms 
are dark, (c) The main tunneling channel from Cu dx2_y2 electrons on nearest neighbor Cu atoms to the tip through overlap 
with the central Cu atom, followed by hopping to the pz orbitals of the apical 0-atom and the Bi atoms, (d) Computed and 
measured spectra in the low energy region marked by gray shading in (a). Inset gives the real (solid lines) and imaginary 
(dashed lines) parts of the self-energy applied to the Cu d^2_y2 states and to all other states. 



where the density matrix pss' = -^^a^ta^'' cxG~^, 
is, in fact, the spectral function written in terms of 
retarded/advanced Green function and the self-energy. 
Eq. [3] differs from the more commonly used Tersoff- 
Hamann approach [l6| in that it takes into account the 
details of the symmetry of the tip orbitals and their over- 
lap with the surface orbitals. 

The use of the spectral function recasts Eq. [3] into the 
form 

^ = E^*"' (4) 

where 

= - f E PtAEF)Vt>sGt^E\G-^,V^,„ (5) 

t' ss' 

and the Green functions and the self-energy are evalu- 
ated at energy E = Ep + eV^. Eq. [5] is similar to the 
Landauer-Biittiker formula for tunneling across nanos- 
tructures (see, e.g., Ref. ITtI ). and represents a slight re- 
formulation of Refs. ^andflgl. 

The nature of Eq. [5] can be understood straightfor- 
wardly: Gscx gives the amplitude with which electrons 
residing on the a^^ orbital in the solid propagate to the 
surface at energy E broadened by T," a- The term Vst 
is the overlap between the surface orbital and the tip, 
while ptt' gives the available states at the tip. Hence, 



Tta gives the contribution of the a^^ orbital to the cur- 
rent, and the summation in Eq. [4] collects these individual 
contributions to yield the total tunneling current which 
reaches the tip. Thus, selecting individual terms in Eq.[5] 
provides a transparent scheme to define tunneling paths 
between the sample and the microscope tip. 

Fig. [31(a) shows the tunneling spectra over the broad 
energy range of ± 1 eV. At high positive voltages, the 
computed spectrum (black line) is fairly structureless. 
At low energies, a gap accompanied by the characteris- 
tic peak-dip-hump features is observed. The calculations 
show the anti-bonding (VHA-a) and bonding (VHS-b) 
VHS's[3] as distinct structures, followed by a broad dip 
around -0.7 eV and subsequent rise near -1 eV. In all 
these respects the present calculations follow the experi- 
ment (red line) of Ref. [1]. Moreover, theory reproduces 
the observed asymmetry of the tunneling spectrum with 
excess intensity at negative biases. The rapid increase in 
current at high binding energies results from increasing 
spectral weight of Cu and other orbitals contributing 
to the 'spaghetti' of bands starting around 1 eV binding 
energy (see Fig. 2(a)). We emphasize that the LDOS 
of the Cu dx'2-y2 (green line in Fig. 3(a)) does not pro- 
vide a good description of the spectrum. In particular, 
the Cu dx2_y2 LDOS possesses an asymmetry which is 
opposite to that of the tunneling spectrum. 

Fig. 3(d) gives a blow up of the low energy region 
of ±0.2ey, shown by gray shading in Fig. 3(a). The 
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computed spectrum is seen to reproduce the coherence 
peaks and the characteristic peak-dip- hump feature. The 
generic form of the real and imaginary parts of the self- 
energies applied to the Cu dj.2_y2 orbitals (solid and 
dashed blue lines, respectively) and the rest of the or- 
bitals are given in the inset. Fig. 3(b) shows the com- 
puted 'topographic map' of the BiO surface in constant 
current mode. Bi atoms appear as bright spots in accord 
with experimental observations, while O atoms sit at the 
centers of dark regions. 
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FIG. 4: Main frame: Partial contributions to the tunneling 
current from various orbitals in the two cuprate layers. The 
Cu02 layer closest to the tip is identified as layer 1 or LI, 
while the second layer is denoted by L2. Specific contributions 
are: (1^2 _y2 orbitals of the four nearest neighbor Cu atoms 
(Ll-nn, blue line); (1^2 _y2 orbitals of the four next nearest 
neighbor Cu atoms of the first layer (Ll-nnn, green line) ; 
d^2_y2 orbitals of the Cu atoms of the second layer (L2, red 
line); of the central Cu atom of LI (magenta line). Inset: 
Decomposition of the current from the second cuprate layer: 
Total contribution (red line); contribution of the four nearest 
neighbours (blue line) ; and the next nearest neighbours (green 
line) . 



An analysis of the partial terms of Eq. [5] reveals that 
the dx2_y2 orbital of the Cu atom lying right under 
the Bi atom gives zero contribution to the current [iot . 
The dominant contribution to the spectrum comes from 
the four nearest neighbor (NN) Cu atoms as indicated 
schematically in Fig. 3(c). A detailed decomposition of 
Eq. m is shown in Fig. [H where paths starting from the 
Cu02 layer closest to the tip (LI), as well as from the sec- 
ond cuprate layer (L2) are considered. The signal from 
cuprate layers is dominated by the dx2_y2 orbitals on the 
four nearest neighbour (nn) Cu atoms in LI up to about 
-0.7 eV (blue line). At higher binding energies, the con- 
tribution from the d^^ electrons from the Cu atom right 
below the Bi atom or the tip grows rapidly (magenta 
line). 

A smaller but still significant contribution comes from 
the four next nearest neighbour (nnn) dx2_y2 orbitals in 



LI spread over a wide energy range (green line, main fig- 
ure), while the total current originating from the dx2_y2 
orbitals of L2 is quite localized over zero to -0.6 eV bias 
(red line, main figure). Fig. [H emphasizes the nature of 
the current associated with the cuprate layers and points 
out an intrinsic electron-hole asymmetry originating from 
the d^-^ orbitals. We note however that the Bi and O or- 
bitals in the surface Bi-0 layer can also play a role in 
producing an asymmetric background current. 

In conclusion, we find that STS spectrum for Bi2212 is 
strongly modified from the LDOS of dx2_y2 by the eflFect 
of the tunneling matrix element. Much of the observed 
asymmetry of the spectrum can be explained within the 
conventional picture due to the turning on of Cu d^^ and 
other channels with increasing (negative) bias voltage. 
This indicates that the eflFects of strong electronic corre- 
lations on the tunneling spectrum are more subtle than 
has been thought previously. However, we should note 
that we have not analyzed spectra associated with the 
deeply underdoped regime where charge order has been 
reported [2l|. The present method naturally allows an 
analysis of the tunneling signal in terms of the possible 
tunneling channels and the related selection rules. Our 
scheme can be extended to incorporate eflFects of impuri- 
ties and various nanoscale inhomogeneities by using ap- 
propriately larger basis sets in the computations. 
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